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Abstract 

In order to explain light curve (LC) for Supernova (SN) we derive 
a classical formula for the conversion of the flux of kinetic energy into 
radiation. We then introduce a correction for the absorption adopting an 
optical depth as function of the time. The developed framework allows to 
fit the LC of type la SN 2005cf ( B and V ) and type lip SN 2004A (B,V,I 
and R ). A relativistic formula for the flux of kinetic energy is also derived 
in terms of a Taylor expansion and the application is done to the LC of 
CRB 050814. The decay of the radioactive isotopes as a driver the LC for 
SNs is also reviewed and a new formulation is introduced. The Arnett’s 
formula for bolometric luminosity is corrected for the optical depth and 
applied to SN 2001ay. 

Keywords: supernovae: general, supernovae: individual (SN 2005cf), super¬ 
novae: individual (SN 2004A ), supernovae: individual (SN 2001ay), supernovae: 
individual (SN 1993J ) 


1 Introduction 


The term light curve (LC) for Supernova (SN) usually denotes the behavior of 
the apparent/absolute visual magnitude as function of the time. The devel¬ 
opment of the multiwavelength astronomy fixes the wavelength passband, i.e 
BVRI, or the frequency zq i.e. 15.2 GHz, or the energy, i.e. lkev. Further on in 
gamma, X and radio astronomies the flux or t he coun t rat e are recorded rather 
than the magnitude, see as an example Fong et all ( 2012 1. The first model to 
be considered is connected with the radioactive decay 


L = L a . oexp(-) 

7"n 


( 1 ) 


where L and La.o are the luminosities at time t and at t = 0 respectively, A is the 
considered wavelength and r n is the typical lifetime, see Bowers and Deemin g 
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( 19841 ). On introducing the apparent magnitude m \, the previous formula be¬ 
comes 


t 


m\ = k ^ + 1.0857(—) 


( 2 ) 


where k' x is a constant. 


The most important radioactive isotopes are 56 Ni 
with T n = 8.757 d and 56 Co with r„ = 111.47 d. The analysis of many au¬ 
thors has shown that the decay of one of these two r a dioac tive isotopes fit 
only few days of a typical LC, see Smith and McCravl ( 20071 ). At the same 


time the spectral i ndex in the radio of SN 1993J is constant after 700 days, 
see Figure 8 in Martf-Vidal et all ( 2011 ) and this observational fact points to¬ 
ward the presence of synchrotron emission having flux , F{y) oc v~ a with 
a ss 0.7. The hypothesis of the synchrotron emission in SNs is not n ew and 
we now report some applica ti ons among oth e rs: GR Bs, see lPreece et al. ( 2002 ): 
iBeniamini and Piran (12013 ): Burgep_et_alJ (2014 ) and Supernovae Remnants 
(SNRs). see iKatsuda et a" ( 2010h : iMiceli et al.l ( 20131) . The presence of the 
synchrotron emission makes attractive the analysis of a turbulent cascade from 
the large scale to the small scales where presumably the relativistic electrons are 
accelerated. Insofar we have isolated two completely different physical mecha¬ 
nism for the source of radiation in the LC: (i) the number of radioactive isotopes 
as function of the time, (ii) the flux of mechanical kinetic energy which is the 
driver for the power injected in the turbulent cascade. The fact that the spectral 
index in the optical regime varies considerably with the time points toward a 
variable optical thickness as function of the time and of the considered pass- 
band. The basic idea is that the optical thickness is low at the beginning of the 
LC and it increases its value with time. A series of questions can now be posed 


• Can we build a formula for the flux/magnitude versus time relationship in 
the framework of conversion of the mechanical luminosity into radiation? 

• Can we introduce the correction for variable optical thickness introducing 
a dependence of the optical thickness with time? 

• Can the new developed framework be applied to the various LCs which 
arise from the various typologies of LCs such as type la, lb, Ic or type lib, 
II-L, II-p, Iln? 

• Can the radioactive model in it’s various versions model the most common 
LCs? 


2 Preliminaries 

Here we first introduce an elementary equation of motion and then we assume 
a linear relationship between the mechanical and the observed luminosity at a 
given frequency v. 

2.1 The simplest equation of motion 

The equation of the expansion of a SN in the first ten years can be modeled by 
a power law of the type 

m = Mh a * ( 3 ) 

to 
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where R is the radius of the expansion, t is the time, Rq is the radius at t = to 
and a is an exponent which can be found from a numerical analysis. The velocity 
is 

V(t) = aRo(^) a t^ . (4) 

to 

As an example in the case of SN 1993J we have a = 0.828. The rate of transfer 
of mechanical energy, L m is 

L m(t) = ^p(t)£nR(t) 2 V(t) 3 . (5) 

We now assume that the density in front of the advancing expansion scale as 


p(t)=po(^-r 


( 6 ) 


where d is a parameter which all ows to match the observations; this assumption 
is not new and as an example [Nagy et ahj ((2014) quotes d=3. The mechanical 


luminosity for the power law dependence of the radius becomes 

L m (t) = L 0 (^) 5a - da ~ 3 , 

to 


(7) 


where Lq is the luminosity at t = to. 


2.2 The emitted radiation 


The energy fraction of the mechanical luminosity, L u , deposited in the frequency 
v is assumed to be proportional to the mechanical luminosity through a constant 

€-V 

Lu — 6iyL m . ( 8 ) 

The flux at frequency v and distance D is 


F v 


;L r 


AttD 2 


(9) 


The problem of the absorption can be parametrized introducing a slab of optical 
thickness r nu . The emergent intensity /„ after the entire slab is 

r= r s^dt , (io) 

Jo 

where S u is a uniform source function. The integration gives 


I v = S v (l-e~ T “) , (11) 

see formula 1.30 in Rvbicki and Liglitmanl ( 199lh . A first model for the optical 
thickness assumes a power law dependence 


_ „ j-a2,v 
— Ui i/t 


( 12 ) 


where ai, v and a 2 ,v are two coefficients which can be found from the astronom¬ 
ical data. The flux corrected for absorption in the power law case, F v , is 


v ' e = 4^D2 


(1 — e ~ T ") 


(13) 
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An expression for the flux of the first model can be obtained inserting the sim¬ 
plest equation of motion for the R — t dependence in the mechanical luminosity 
as given by 0 

Fv, c = F Vt o(^-) 5a ~ da ~ 3 (1 — e~ Tu ) , (14) 

to 

where F v fi is the flux at = to- This formula is useful when we have the 
flux versus the time, as an example Jansky versus JD. The absolute/apparent 
magnitude version of for the first model is 


m(t) 


2.5 In (t) ad — 12.5 a In (f) + 7.5 In (t) — 2.5 In (l — e ai *° 2 ) 
In (2) + In (5) 


+ m k , 
(15) 


where m k is a constant of calibration. This formula is useful when we have the 
absolute/apparent magnitude versus time as in the case of optical LC in SN. 
The asymptotic approximation is 


m(t) 


1.085ad In (t) 


5.428a In (t) + 3.257 In (f) + rrik + 


0.542 



• ( 16 ) 


A second model for the optical thickness assumes an exponential law depen¬ 
dence 

r„ = a Uv (l —, (17) 

where a, 2 ,v and 03 ^ are three coefficients which can be found from the 

astronomical data. The absolute/apparent magnitude version for the second 
model is 


2.5 In (t) ad 


12.5a In (t) + 7.5 In (t) — 2.5 In 1 — e 


—ai t, ( 1—e “ 2 ’’ 

— e v 


In (2) + In (5) 



+m k ■ 


(18) 


A third model for the optical thickness makes a comparison between the 
observed intensity I 0 b s and the theoretical intensity Ith through the optical 
depth 

1 _ e - T - = ^ . (19) 

Uh 

The optical depth is 

= — ln(l — —j—) ■ ( 20 ) 

hh 

The observed intensity as function of the time is an astronomical quantity and 
the theoretical intensity can be the mechanical luminosity or the momentary 
number of radioactive isotopes. Once the temporal behavior of r v is derived we 
search for the best fit as function of time. A fit already used is the power law fit 
as represented by eqn. m- Another type of fit is the logarithmic polynomial 
approximation of degree M, 

T„(t) = a 0 + oi(ln(i)) + a 2 (ln(t )) 2 J-1- a M (hi {t)) M . (21) 

The presence of the logarithm allows to cover the oscillatory behavior of r„ over 
many decades in time. 
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Figure 1: The V LC of SN 2005cf (empty stars) and theoretical curve as given 
by the first model, see eqn. (fl5l) (full line). 



Figure 2: The V LC of SN 2005cf (empty stars) and theoretical curve as given 
by the asymptotic formula (flTH) . (full line). 


3 Astrophysical results 


The time is usually expressed in JD or seconds and a subtraction of the initial JD 
or seconds relative to the considered phenomena should be done in order to have 
zero at the beginning of the temporal scale. We start by analyzing SN 2005cf 
in NGC 5812 which is of type la, it’s distan ce is 29.4 Mpc and the distance 

Figure U] reports the temporal 


moduli n = 32.51, see (Pastorello et al 


evolution of the V visual magnitude of SN 2005cf for the power law model as 
well the interpolating curve and Figure [2] the asymptotic approximation; data 
as in Table [I] 

The quality of the fits is measured by the merit function y 2 


X' 


— 'y v (mth nt obs ) 


where mth and m 0 b s are the theoretical and observed magnitude, respectively. 

The (B-V) color evolution of SN 2005cf for the exponential law model (data 
as in Table [2]) is reported in Figure [3] 

The SNs of type Up are characterized by a flat LC for a long period of 
time, i.e. 100 days. We therefore analyzed SN 2004A in NGC 6207 , which is 
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Table 1: Numerical values of the adopted parameters and \ 2 for the optical LC 
in SN in the case of optical thickness with a power law dependence, a = 0.828 
everywhere. 


Name SN 

band 

d 

ai 

a 2 

m k 

x 2 

SN 2005cf 

V 

3.80 

1.9 10" 4 

2.95 

6.73 

0.188 

SN 2005cf 

B 

3.85 

1.98 10" 4 

3.2 

7.11 

6.7 

SN 2004A 

V 

4.15 

1.0 10" 4 

4.15 

5.33 

12.4 

SN 2004A 

B 

3.52 

2.6 10~ 3 

1.85 

9.53 

2.50 

SN 2004A 

I 

4.33 

2.0 10~ 5 

2.65 

3.4 

0.35 

SN 2004A 

R 

4.05 

7.8 10“ 5 

2.3 

4.81 

1.10 


Table 2: Numerical values of the adopted parameters and 'x for the optical LC 
in the case of optical thickness with an exponential law dependence, a = 0.828 
everywhere. 


Name SN 

band 

d 

ai 

a 2 

03 

m k 

x 2 

SN 2005cf 

V 

3.79 

10 

1.75 10 _b 

3 

6.73 

0.202 

SN 2005cf 

B 

3.81 

5 

8.2 10~ 5 

3 

7.32 

6.45 

SN 2004A 

V 

4.36 

10 

6.27 10 -7 

3 

4.07 

1.70 

SN 2004A 

B 

3.79 

100 

6 10" 7 

3 

8.15 

2.4 

SN 2004A 

I 

4.67 

2.5 

1.0 10" 6 

3.0 

1.62 

0.52 

SN 2004A 

R 

4.23 

10 

4.51 10 -7 

3 

3.83 

1.6 


of type Up, the dis t ance is 25.6 Mpc and the distance moduli /z = 31.99, see 
Hendrv et, al. 1 2006 1: iTsvetkov (2 0081 ). Figure 0] reports the temporal evolution 


of the V visual magnitude of SN 2004A for the exponential law model as well 
the interpolating curve, data as in Table [2] 

We now apply the devel oped theory to m o del th e radio flux den sity of 


SN1993J at 15.2 GHz, 
data available at 


see 


Poolev and Green (1993); Ho et al 


with 


http://www.mrao.cam.ac.uk/~dag/snl993j.html . In this radio-case we plot 
the flux version of the first model as given by eqn. Q. see Figure [5] and Table 
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Table 3: Numerical values of the adopted parameters and y 2 for the radio LC 
in SN 1993J and gamma LC in GRB 050814in the case of optical thickness with 
a power law dependence, a = 0.828 everywhere. 


Name 

band 

d 

Ol 

02 

Fu,o 

x 2 

SN1993J 

15.2 GHz 

2.22 

8.5 10" 4 

1.86 

1.64 

7.9 10” 4 

GRB 050814 

0.2-10 kev 

2.79 

0.026 

1.259 

8.3810 -08 

4.9 10" 18 


The theory is now appli ed to G RB 050814at 0.3-10 kev in the time interval 
10 -5 — 3 days, see Jakobsson et al. (,20061 ) with data available at 


6 


























































Figure 3: The (B-V) color evolution of SN 2005cf (empty stars) and the relative 
fitting straight line (full line). The theoretical curve is given by the second 
model. The time is limited to the first 25 days. 



Figure 4: The V LC of SN 2004A (empty stars) and theoretical curve as given 
by the second model, see eqn. (fl9l) (full line). 



Figure 5: The radio flux density of SN 1993J at 15.2 GHz (empty stars) and 
theoretical curve as given by the first model, see eqn. (fT4l) (full line). 
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Figure 6 : The XRT flux of GRB 050814at 0.2-10 kev (empty stars) and theo¬ 
retical curve as given by the first model, see eqn. (1141) (full line). 


http://www.swift.ac.uk/xrt_curves/00150314/. Figure [ 6 ] reports the LC, in this 
case the flux, as function of the elapsed time since Burst Alert Telescope (BAT) 
trigger and Table [3] reports the involved parameters. 


4 Relativistic model 

The density, p, of the ISM at a distance r from the SN is here modeled by a 
Lane-Emden (n = 5) profile 

*> = *(!+ £)./» ’ <22) 

where b represents the scale. The relativistic conservation of momentum for the 
thin layer approximation in presence of a the Lane-Emden (n = 5) profile is 
given by the following differential equation 

4=b 3 {r(t)) 3 pnV3f t r(t) _ Ab 3 r 0 3 pn V3/3 0 

(3 b 2 + (r(t)) 2 ) 3 / 2 c\J —isTWll + i (3 b 2 + r 0 2 ) 3 / : V-A)“ + 1 

where ro is the initial radius of the advancing sphere, vq is the initial velocity 
at r 0 , c is the light velocity and 0o = qr- The relativistic transfer of energy 
through a surface, A, is 


L m , r = Ay 2 (pc 2 +p)i 


(24) 


where p is the pressure, for sake of simplicity we take p=0, and the Lorentz 
factor 7 is 

7 "" v'l - P 2 dt ’ ^ 

see eqn. A31 in De Young! ( 2002 ) or eon. 143.441 in Mihalas and Mihalas ( 2013 b 
In the case of a spherical cold expansion 


L m ,r = 4t rr(f ) 2 


-p(t)c 3 f3(t) 


1 - /3(t) 2 


(26) 



















Figure 7: Numerical L m ^ r computed according to Eq. (051) (full line) and series 
solution as given by Eq. (051) (dotted line). Data as in Table [I] 


We now assume the following power law behavior for the density in the advanc¬ 
ing thin layer 

p(t) = Po( t j) d , (27) 

and we obtain 

Lm,r = 4tt r{tf 1 _ * 2 p 0 (j) d c 3 P{t) . (28) 

We can now derive L m ^ r in two ways: (i) from a numerical evaluation of r(t) 
and v(t), (ii) from a Taylor series of L mjr (t) of the type 


3 

L m ,r{t ) = ^2 Gn ~ ^ 
71=0 


(29) 


The coefficients are 


a 0 = 


ai = 


a 2 = 


«3 = 


-4Trr 0 2 (-^) d c 3 /3 0 

4 7rr 0 c 4 ffo 2 (^) d (9b 2 /3o 2 +3fe 2 -2r 0 2 ) 

(3 b 2 +r 0 2 )(/3o 2 -l) 

2 /3 0 3 c 5 (^) d 77 (l62 b 4 /3 0 4 —297 b 4 /3 0 2 -9 b 2 /3 0 2 r 0 2 -45 b 4 +15 b 2 r 0 2 -2 r 0 4 ) 

(3 b 2 +r 0 2 ) 2 (/3o 2 —1) 

18 (270 b 2 /3o 6 -675b 2 / 3 o 4 -33/3o 4 r o 2 +480 b 2 /3 0 2 +62/3 0 2 r 0 2 +45 b 2 -5 r 0 2 )/3 0 4 b 4 c 6 7i- (^) d 
r- 0 (27b 6 +27b 4 ro 2 +9b 2 r 0 4 +ro 6 )(/3o 2 -l) 


(30) 


Figure [7] compares the numerical solution for the luminosity and the series ex¬ 
pansion for the luminosity about the ordinary point to- 
The flux at frequency v and distance D is 




47T.D 2 


(31) 


The flux corrected for absorption in the relativistic case is 


F — 

^ i',c,r — 




47tD 2 


(1 ~ e~ Tu ) 


(32) 
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Table 4: Numerical values of the parameters used in relativistic solutions. 


parameters 

t 0 = 2 -8 yr or to = 0.63 s ; ro = 0.00195 pc ; /3 q = 0.833 ; b = 0.004 pc 



Figure 8: The XRT flux of GRB 050814at 0.2-10 kev (empty stars) and the¬ 
oretical curve as given by the relativistic numerical model, see eqn. (1321) (full 
line). 


As a behavior for t v as function of time we select a logarithmic polynomial 
approximation, see (ecinl2Tl). of degree 9 and Figure [5] reports the flux of the 
relativistic LC as function of the elapsed time since Burst Alert Telescope (BAT) 
trigger. 


5 The radioactive model 

Here we consider the decay of a radioactive isotope, a radioactive chain and the 
Arnett’s rule for the bolometric luminosity. 


5.1 Decay of one element 

The decay of a radioactive isotope is modeled by the following equation 

-dN = —dt , (33) 


where r n is a constant and th e negative sig n ind icates that dN is a reduction in 


the number of nuclei , see Yang and Hamilton ( 2010l) . The integration of this 


differential equation of the first order in which the variables can be separated 
gives: 

N(t) = N 0 e~^ , (34) 
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where Nq is the number of nuclei at t = 0. The half life is T) / 2 = ln( 2) r n . The 
absolute magnitude version of the previous formula is 

M = -C Log w {N{t)) =—- + k , (35) 

Tn 

where M is the absolute luminosity, C and k are two constants. This means 
that we are waiting for a straight line for the absolute magnitude versus time 
relationship. At the same time the observational fact that the spectral index 
in the radio varies consid erab l y bu t becomes constant, /3 « —0.7, after ss 700 
days, see Figure 8 in Marti-Vidal et ajJ ([20111 ). asks the absorption. 


5.2 Radioactive chains 

The isotope 56 Ni is unstable and decays (n = 8.757 d, 7 \/ 2 =6.07 d) into 56 Co 
emitting gamma photons. The isotope 56 Co is unstable and decays (r 2 = 111.47 
d, T x / 2 = 77.27 d) into 56 Fe through electron capture and /3-decay. The decay 
rates of of the two species, species 1 is 56 Ni and species 2 is 56 Co, is modeled 
by the following equations 


djv 2 (t ) = M _ Ml. 

df ri r 2 


(36a) 

(36b) 


The two solutions obtained inserting as initial conditions iVi(0) = and 
iV 2 (0) = 0 are 


Ni = iVope r i 

( T 2 Nq,1 _L 


N 2 = 


^2^0,1 


1 \ — 


V T 1 ~ T2 


Tl - T 2 


(37a) 

(37b) 


The_sum of the two species, N(t), is according to formula (8.5) in IRust et al.1 

(l2010h 

N(t)=C 1 N 1 (t) + C 2 N 2 (t) , (38) 

where C\ and C 2 are two adjustable parameters. This linear sum is associated 
with the LC in SNs assuming that the 7 -rays are thermalized in the ejecta and 
emerge in the various bands. The logarithmic form,M(t) , is associated with the 
magnitude evolution 

M = ^ (39) 


MD 


where 


— ln(- 


1 


MN = fcln(2) + fcln(5) — ln(JVo,i) 

_ t_ t(rx — r 2 ) _ t_ _ t_ _ t_ 

-(C 2 r 2 e T 2 C u -2 + Cie un-Cie n t 2 - C 2 r 2 e ^)) , (40) 


ti - t 2 
where k is a constant and 


Md= ln(2) + ln(5) 


(41) 


We plot the decay of the LC of SN2001el , which is of type la, adopting a 
distance modulus of 31.65 mag, see iKrisciunas et ah (2003). the nuclear decay 
which according to equation (1551) is a straight line, and the theoretical curve of 
the two species as represented by equation l39l see Figured] 
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Figure 9: The V LC of SN 2001el (empty stars) in absolute magnitude, the 
theoretical curve as given by equation (1351) when the radioactive decay of the 
isotope 56 Ni (r„ = 8.757 d or Ti/ 2 =6.07 d , fc=-18.65) was considered (full 
line), and the theoretical curve, (dashed line), connected with the decay of 
two species as represented by eg. (1391) (n = 8.767d, r 2 = 111.477d, iVop = 1, 
( 71 = 1 . 9995 , <7 2 =1.0005 and fc=-18.4). 


5.3 Bolometric luminosity 

The bolometric LC after I Arnett (' 1982 1 : lArnett et~ak (1985) has been associated 


with the combined radioactive decays of the isotopes 56 Ni and 56 Co. A formula 
of practical use is given by 

L(t R ) = a ((6.45 x 10 43 )e ~ tB/8 ' 8 + (1.45 x 10 43 )e“ tB/111 - 3 ) — , (42) 

where t R is the elapsed tim e from the ex plosion to the maximum of the LC and 
a is ~ 1, see formula (2) in Krisciunas et all ( 2011 ). The previous formula rep¬ 
resents the optically thin case. According to the comparison method developed 
in Section [2.21 a corrected bolometric luminosity for the absorption, L c (t R ), is 


L c (t R ) = L(t R ) x (l —e- ait “ 2 ) , (43) 

and Figure [TO] reports the comparison between observed and theoretical bolo¬ 
metric luminosity. 


6 Conclusions 

Classical and relativistic flux of energy: The classical flux of kinetic energy can 
be easily parametrized in the case of a radius-tinre relationship represented by a 
power law, see eqn. 0. Conversely is more complex to derive the relativistic flux 
of kinetic energy which requires a relativistic law of motion. In the framework of 
Lane-Emden (n = 5) profile as given by ean. (l22l) and momentum conservation 
in a thin layer we can deduce an analytical solution for the relativistic flux of 
energy in terms of a Taylor series, see the four coefficients in (1301) . 

Light curve: Assuming a linear relationship between the luminosity in the 
various astronomical bands and the classical or relativistic flux of mechanical 
kinetic energy we can easily deduce a theoretical time dependence for the LC, 
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Figure 10: The bolometric LC of SN 2001ay (empty stars) and the theoretical 
curve as given by eqn. (full line). The parameters are tn=22 days , a = 1 
and the comparison method with a power law fit gives aq = 

— 1.33 5. The astronomical data are extracted from Figure 18 in 
(201 ll ) by the author. 


105.88 arid q 2 = 
Krisciunas et al. 


see classical eqn.(0 or relativistic ecm. . This theoretical dependence is not 
enough and the concept of optical depth should be introduced. Among the 
infinite relationships for optical depth as function of time we selected a power law 
dependence, see eqn. m , an exponential behavior, see eqn. ED, or a logarithmic 
polynomial approximation, see eqn. ED- 

Nuclear Decay: The LC of a SN is often model by the decay of the radioactive 
isotope 56 Ni, but in order to follow the LC with time we should introduce a 
radioactive chain, see eqn. Further on some classical approach to the 

bolometric luminosity must be corrected for the optical depth, see ean. (|43D . 

Comparison with astronomical data: The framework of conversion of the 
classical flux of mechanical kinetic energy into the various optical bands coupled 
with a time dependence for the optical depth allows to simulate the various 
morphologies of the LC: for a type la we chosen SN 2005cf, see Figure [T] for V 
band and Figure [3] for (B-V) color. The enigmatic behavior of type Up SNs, 
here represented by SN 2004A, can also be modeled, see Figure[4]for the V band 
and Table [T] for B, I and R bands. The opposite sides of the electro-magnetic 
spectrum can also be simulated: for the radio band of SN 1993J see Figure [5] 
and for the gamma/X spectrum of GRB 050814 see Figure [G] More complex 
is the derivation of the relativistic flux of energy here parametrized by a series 
expansion. The coupling of the previous series with a logarithmic polynomial 
approximation allows to model fine details such as the oscillation in LC visible 
at ~ 1000 s in GRB 050814, see Figured All the fits here presented report the 
X 2 , see Tables [T| and [2] This means that other types of functions for the optical 
depth versus time have a reference for comparison. 
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